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ABSTRACT 


This is the final report to NASA Contract NAS 6-2571. 
Sources of geometric and radiometric fidelity errors in 
AN/APQ-102A radar imagery are discussed, along with a 
digital computer program to correct the distortions. The 
major effort, a computer program which will process 
digitalized recorded AN/APQ-102A phase histories into 
imagery, is described. All computer programs are listed 
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SECTION I - INTRODUCTION AND SUMMARY 

This report describes the work accomplished on a study program entitled, "Study of 
Synthetic Aperture Radar Imagery Characteristics," funded under NASA Contract NAS 6-2571. 
The objective of the program was to analyze the characteristics of synthetic aperture radar 
(SAR) imagery and develop digital processing techniques to utilize this type of data in con- 
junction with other sensor data on the NASA earth resources program. Specifically, the 
major effort was directed toward producing the digital computer programs for the processing 
of data obtained from the AN/APQ-102A radar system. 

The study program consisted of two major tasks: (1) the definition of SAR image character- 
istics, and (2) the development of digital computer programs to accept phase history data and 
generate a radar image normalized relative to both intensity and geometry. 

Section II discusses the sources and magnitude of errors in the AN/APQ-102A imagery. The 
theoretical analysis consists of enumerating the known sources contributing to geometric 
distortion, determining the effect of each source, and combining to yield an overall estimate 
of the geometric fidelity of the imagery. Sources for geometric distortion fall into three 
categories: (1) sensing geometry, (2) radar equipment errors, and (3) errors in the aircraft 
inertial navigation system (INS) and altimeter. In addition to the theoretical analysis, the 
distortions in an AN/APQ-102A image of Wallops Island flown on 30 August 1973 were 
measured. 

Section III describes the computer programs and procedures developed to process AN/APQ- 
102A phase history data. These programs and procedures were validated by actually 
processing imagery at the Wallops Station facility utilizing the Optronics Microdensitometer 
and the Honeywell 625 computer. This validation effort included the training of Wallops 
Station personnel, thus giving NASA the capability to process subsequent radar data without 
contractor support. 

Conclusions are given in Section IV, and the appendixes contain program listings of all 
computer programs generated or used. 


- 1 - 
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SECTION II - DEFINITION OF IMAGERY CHARACTERISTICS 


1. GENERAL 

The AN/APQ-102A has been quite successful in mapping for tactical purposes; however, 
its imagery has small geometric and radiometric fidelity errors which it would be 
desirable to remove when it is being used for cataloging earth resources. Some of the 
geometric errors are internally generated within the radar; however, these errors are 
generally small. The major sources of geometric errors are inertial system errors. 
Since these errors are not known for any particular flight, their effect (geometric dis- 
tortion) must be measured by comparison with a map or other well -controlled data. This 
section discusses the error sources, their effect on geometric fidelity, and a method of 
measuring geometric distortions through the use of terrain features recognizable both in 
the radar image and on a map. 

The basic design of the AN/APQ-102A includes features that minimize radiometric dis- 

2 l/2 

tortions that would be caused by sensing geometry (e.g. , esc cos vertical antenna 
pattern). Radiometric distortions can be determined by measuring the deviation of the 
radar transfer function from the ideal or by imaging a calibrated radiometric range. 

Only the first of these methods is discussed. 


2. GEOMETRIC FIDELITY ANALYSIS 

The velocity and flight characteristics of the RF-4 aircraft and its avionics are used in 
the numerical calculation of the magnitude of the error components. The calculation is 
typical of the parameters of the flight of 30 August 1973. This analysis includes only 
fixed target imagery and only the modes listed in Table I. 


Preceding page blank nqx 
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TABLE I -AN/APQ-102A HIGH-RESOLUTION OPERATING MODES | 


Mode 

Altitude 

(ft) 

Type of 
coverage 

Range coverage 
(NMI) 

Velocity 

(ft/s) 

1 

500 to 5, 000 

HR 

0 to 10 both sides 

700 to 1250 

5 

30,000 to 50,000 

HR 

5 to 15 both sides 

700 to 2000 

6 

30, 000 to 50,000 

HR 

10 to 30 left 

700 to 2000 

7 

30,000 to 50,000 

HR 

10 to 30 right 

700 to 2000 


3. ACROSS-TRACK ERRORS 


a. Ground Range Sweeps 

The CRT recorder in the AN/APQ-102A employs ground range sweeps. Two 
characteristics of the sweep are normally considered relative to geometric fidelity, 
i.e. , linearity and stability. Sweep linearity is expressed in terms of percent error 
in the distance between two points in the sweep interval. The linearity of the sweep 
of the AN/APQ-102A is ±0.5 percent. The error in position of a target at one edge 
of the swath with respect to the other for a 10-NMI swath is 

±(0.005) X (5) (6080) = ±152 ft . (1) 

The long-term stability of this error should be good, and the error should be highly 
correlated between data films using the same recorder on successive missions. 

The expected dominant spatial frequency of this range scale factor error is one-half 
cycle per sweep length, with the error near the center of the sweep trace being very 
small. 

b. Film Thickness Variations 

Linear film thickness variations can cause errors on the order of four feet in the 
range direction. The spatial frequencies of these errors have not been determined. 

-4- 
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c. Range Displacement Error from Target Altitude 

Conversion from slant range measurement to ground range measurement depends on 
the relative altitude between the radar flightpath and a given target and thus is 
affected by terrain roughness, earth curvature, etc. No error is attributed to the 
radar for this operation. 


4. ALONG-TRACK ERRORS 

a. Recorder Film Drive Error 


The major error in film drive is caused by the error in the measurement of ground - 
speed. The accuracy of the velocity measuring equipment is about nine feet per 
second. Thus, the linear error in target location resulting from errors in velocity 
is 


where 




» 


V_ = error in velocity 
E 

V = aircraft velocity 

D = the along-track distance over which error is to be considered. 
In a high-speed mode over a five-mile distance, the linear error is 


( 2 ) 


9 

1500 


X 5 X 6080 = 180 ft 


(3) 


It can be further assumed that a three-mil peak sinusoidal error is present because 
of eccentricity of the film metering drum. This would be equivalent to a three-foot 
error. 
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b. Film Thickness Variations 

The errors from film thickness variations are estimated to be small (about four feet). 


c. Clutterlock Stability 

The stability of the clutterlock motion compensation loop is about five Hertz and 
would produce an error in the longitudinal (y) direction according to the following 
relationship: 


dy n 



df A 
d 

2 V 

y 



where 

frequency error = df^ - 5 Hz 
wavelength = A - 0. 1 ft 
groundspeed = V = 1500 ft/s 

y 

e = relative error 


( 4 ) 


e . dyj - dy 2 

df A 

e ■= (R . - R ) 

mm max 2V 

y 


30,400 X 5 X 0. I _ 

f - — " 5 1 1 

2 X 1500 

This is a skew-type error, and its frequency is estimated to be very low. 


d. Correlator Film Drive Error for Optically Processed Imagery 

Aspect ratio error contributes a small, steady-state error in the image azimuth 
scale factor. This error is held to less than a resolvable element and is estimated 
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to be 15 feet. Data film drive error resulting from metering drum eccentricity is 
the same as that of the recorder and is three feet. The image film metering drum 
eccentricity permits a maximum sinusoidal error of 23 feet peak at a period of 
3.9 NMI. 


e. Residual Error from Motion Compensation Instrumentation 

Compensation for sideways motion in the AN/APQ-102A radar is achieved in the 
following manner. An antenna -mounted accelerometer system measures the side- 
ways accelerations which are integrated and combined with the clutterlock- 
measured velocities. The clutterlock takes an average of range samples at intervals 
along the entire 10-mile swath. The motion compensation signal thus derived is 
applied at midrange. Thus, no along-track error exists at this point; however, an 
error does exist on each side of the midpoint, with maximum error at maximum 
and minimum ranges. 


To obtain an expression for this error, consider a velocity in the X direction, V . 

X 

A correction is made so that no error exists at Y . However, there is a difference 

o 

between the hyperbola on which mapping is occurring and the straightline correction 
that is applied. From Figure 1, the following expressions may be written : 


Y = (R 2 

0 


2 1/2 

h ) cos 0 


Y = 
o 



R 

o 


( 5 ) 

(6) 



V l/2 

r~ R = (R 2 + h 2 ) cos 0 
V y 0. '0 


( 7 ) 

(8) 
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COS Ip = 


< R o * h2) 


2 2 1/2 V Y ° 

(R + h 2 ) : 

2 2 
(R q + h ) 


2 ,2 1//2 V 

(R + IQ n _2L 

2 2 1/2 ° V Y 

(R + h 2 ) Y 

o 


The along -track error is 


Y - Y = 


Xx Xx R P 2 + h 2 > ' 

V Y V Y ° 2 2 1/2 

Y Y (R + h 2 ) 

o 


= I R - R 


2 2 1 ^ 2 1 V 

(R + h ) I X 


0 2 2 1/2 V Y 

<R 2 + h 2 ) J Y 


! 


(9) 


( 10 ) 

/: 


n 


(ll) 
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If the velocity, V , is 10 ft/s, V is 1500 ft/s, R is 20 miles, R is 15 miles, and 

2\. Y O 

h is 40, 000 ft, the error is 

Y - Y ,= 26.7 ft . (12) 

0 

_f. Effect of Clutterlock and Across -Track Velocity Measurement Error 

An across-track velocity measurement error introduces a squint or skew into the 
final image. An across-track velocity measurement error of 6 ft/s at 1500 ft/s 
produces a pointing error of 

ANGULAR « CTV = = 4 X 10 3 radian . (13) 

The linear error along track of a target on one edge of the swath with respect to the 
other is 


e CTV 


v„ 


(R max “ r min ) 


6 

1500 


(10) (6080) 


- 243 ft 


(14) 


g;. Effect of Vertical Velocity Measurement Error 


An expression relating vertical velocity to along-trnck error can be developed 
similar to that for across-track velocity: 


'vv = h 


<R 2 ± h 2 ) 


1/2 


lK + ” 2 > 


1/2 


(15) 
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For the conditions 


V = 3.5 ft/s 

Lj 

V y = 1500 ft/s 

h = 40,000 ft 

R = 20 NMI 

R = 15 NMI 
o 

therefore, 


e 


VV 


26.9 ft 


h . Effect of Vertical Velocity Measurement Error 

The along -track effect of vertical velocity measurement error may be determined 
from the expression 

£AV z * x AV z • < 16 > 

For the parameters 
AV Z = 2 ft/s 
V = 1500 ft/s 
h - 40, 000 ft 

the along -track error resulting from vertical velocity measurement error is 

40,000 

eAV z = T50T X2 = 53 - 3ft • < 17 ) 
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i. Antenna Pitch and Yaw Errors 


Errors in antenna pitch and yaw will affect along-track geometric fidelity. The 
result is that the clutterlock attempts to correct for the error, causing a skew in 
the imagery. Consider the geometry of Figure 2. The aircraft is flying at velocity 
with an antenna pitch of 8 I and yaw of e & . The error in the along-track direction is 
the mismatch between the best-fit doppler line and the antenna pattern intersecting 
the ground. This error is Y - Y^. 

From Figure 2, the following expressions may be written; 


Y 

ip 



Y = 

Y = 
o 


2 1/2 

-t- h ) cos ib 


Re + he 

a p 


Re + he 

o a p 


(18) 

( 19 ) 

( 20 ) 
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Figure 2 - Geometry of Antenna Pitch and Yaw Error 
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R 0 
o a 


i h 0 ■- 
P 

COS Ip = 


I /u 


(ir i m 

o 


COS i / 1 


R 6 + hd 

o a p 



( 21 ) 

( 22 ) 


The error in along -track position is 


2 2 1/2 

Y - Y = Re + he - (R +■ h ) cos tp 
ip a p 


(23) 


and, substituting for cos ip, the following is obtained: 

, 1/2 


Y - Y. = R0 > he - 
4> a p 


2 2 

R ^ h 

2 2 

R + h 
o 


(R e *■ he ) 
o a p 


(24) 


The error at R , Y is zero but increases as R increases or decreases from R . 

o o o 

Therefore, R q is placed at the midswath position, and R may vary five NMI on either 
side. 

For the following conditions, 


R = 20 NMI 

R = 15 NMI 
o 

h = 40,000 ft 
e = 0. 25 deg 

3 

0 =0.25 deg , 

P 

an along -track error of 90 ft occurs. 
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j. Aircraft Turning Error 

The inability of the aircraft to fly a perfectly straight path introduces errors in the 
along-track direction. When the aircraft goes into a slow turn within the time con- 
stants of the clutterlock, the clutterlock is able to keep the physical beam oriented 
on the zero doppler line. However, the radar thinks it is following a straight line 
flightpath and, as a result, the imagery is skewed. The geometry of this situation 
is depicted in Figure 3. Here, r is the uncompensated turning radius of the aircraft, 
R is the range of interest of the radar, s is the length of the ground track in rotating 
through the angle 0, and S is the distorted, recorded length of the flightpath. The 
error caused by the aircraft turn may be written as 


S - s _ (r f- R)0 - rfl R 
T s r0 r 


(25) 



4779-4 


Figure 3 - Geometry of Errors Resulting from Uncompensated Turn 
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It is difficult to determine what the minimum uncompensated turning radius for the 
R F-4 or C-54 is, but the manufacturer of the UK -4 has indicated that it will fly 
straight with a minimum turn radius of 2050 NMI. Therefore, the error resulting 
from a turn expressed in percent is 

T> IQ 

<r T = — = ~ - - 0.00378 - 0.38 percent . (26) 

When one considers a strip five NMI long, the turning error is 
e T = (0.0038) X 5 X 6080 = 115 ft . 

5. MEASUREMENT OF GEOMETRIC ERRORS 

a. General 

It can be seen from the foregoing that most of the fidelity errors are random and 
are irreducible prior to imaging, because they represent the utilization of onboard 
sensors and their attendant errors. Certain errors such as residual motion compen- 
sation error and the skew introduced by use of an offset frequency for demodulation 
could be removed after flight if the offset frequency and the three-axis translation 
of the aircraft were recorded during flight. These data, however, are not available 
in the AN/APQ-102A, and hence the aforementioned errors are not separable from 
the random errors. 

As a postflight procedure, the processed radar image may be rectified by use of a 
computer program which does a least-square-error fit using precisely imaged 
points whose geographic corrdinates are accurately known. Care must be taken that 
points utilized are coincident with the points whose coordinates are known. 


- 14 - 
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b. Program Rat-ionnlo 

For the distribution correction program, the following radar image distortions were 
considered: 

1. Scale (range versus track) - caused by the separate scaling 
mechanisms involved (in the radar) in the track and range 
directions 

2. Skew - caused by either radar antenna or correlator slit 
misalignment (results in range and track nonorthogonality) 

3. Residual distortion - that which remains after items 1. and 2. 
have been accounted for (caused by height differences between 
the radar ground plane and the elevation of objects in the ground 
area being imaged, nonlinearities in the recording CRT, 
measurement errors, etc,). 

c. Program Steps 

Since the program determines the d istortion relative to (radar) range and track, 
a prerequisite for the analysis is that the radar image measurements (of identified 
ground control points) be performed with the X-axis of the measurement device 
aligned with the track direction of the radar. The steps taken by the program in 
analyzing distortion are summarized as follows: 

1, The ground coordinates of the control points are preliminarily 
aligned with the image coordinate system. This is done by 
determining the relative orientation of two designated control 
points in both the image and ground frames and rotating the 
ground system to coincidence. For correctly signed printout 
of scale, skew, and distortion, it is desirable to choose the 
ground axes to lie within 45 degrees of the image coordinate 
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system (for an existing system, this is accomplished by 1 

controlling the sign and (X-Y) designation of the axes). The 

designation of alignment points (as a program parameter) 

prevents the use of points known to have a high probability of 

substantial relative distortion (e.g. , points extremely close 

together). The selection of two widely spaced points will 

I ' 

suffice for the initial alignment 

2. A least-squares fit between image and ground range coord i- 
ates is performed. The range errors remaining after the fit 
are computed, and the linear correlation coeffficient between 

range errors and track coordinates is determined (if a nonzero p 

coefficient exists, it indicates a residual misalignment). The 

ground coordinates are then rotated to make the coefficient 

zero. This process prevents individual control point errors 

from introducing substantial alignment errors 

r 

3. A least-squares fit between the range image and ground £ 

coordinates is computed and residual errors determined 

(at ground scale) 

4. Track image and ground coordinates are scaled via a least- 
squares fit and average error determined. Skew is then 
introduced into the image coordinate system via two equations: 

Y = Y (27) 

and 

x’(I) = X(I) + A • Y(I) , (28) 

U 
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where 


(X, Y) = original image coordinates 

r t 

(X ,Y ) = skewed coordinates 
A = tangent of skew angle. 

The skew angle is varied (in sign and magnitude) until track 
errors are minimized (as measured by successive least- 
squares fits) 

Several types of analyses are then performed by the program 
to demonstrate the relative contribution of various error 
sources. In each of them, the residual error variance and 
individual point errors (ground scale) are computed (and 
displayed for examination) after various types of image 
correction are introduced. The four types of correction are: 

a. A magnification equal to the average (range and 
track) scaling difference between image and 
ground coordinates 

b. Differential scale correction 

e. Magnificat ion plus skew correction 

d. Differential scale and skew correction. 

The image range scale, track scale, and skew are displayed. 
Plots are created (using CalComp software) which illustrate 
the ground position and residual error of the control points 
after the various types of correction. 
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d. Program Parameters 

The program listing is given in Appendix A. 'ITie utilization of the results to 
restitute digitally processed imagery is discussed in Section III. 

Program parameters are as follows: 

ALPHA 1 - ALPHA 4 Analysis section headers 

N Number of control points 

Nl, N2 Number (corresponding to program 

order) of point pair to be used for 
initial alignment. Nl and N2 should 
be ordered so that Nl has a smaller 
track coordinate than N2 (however, 

Nl < N2 need not be true; i_.e. , 
control points may be entered in 
any order) 

SCPLT Scale of plots (ground units/inch) 

which are created of control point 
positions 

SCERR Scale of error vectors 

(P(I), Q (I)) Ground coordinates of control point 

(X(I), Y (I)) Image coordinates of Ith control point 

A (I) Alphanumeric control point 

designator. 


B 


n 


11 


- 18 - 
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e, Results 

A section of the AN/APQ-102 data film which was flown for NASA on the 30 August 
1973 flight was optically correlated on a laboratory correlator, and the portion of 
the imagery around Wallops Island was analyzed. Sixteen points were used, of 
which coordinates for 6 were from triangulation and 10 from a map. After distor- 
tion correction, the mean error along track was zero with a standard deviation of 
77.52 meters and mean across-track error of 7.2 meters, with a standard 
deviation of 87.7 meters. 

6. RADIOMETRIC ERRORS 

a. General 

As mentioned previously, the transfer function of the AN/APQ-102A was designed 
to compensate for sensing geometry, ideally resulting in no radiometric distortion. 

To the extent that the radar transmitted power and the receiver gain remain constant 
(this can be accomplished by disabling the automatic gain control), the radar can be 
designed to compensate for changes in slant range and depression angle. These 
compensations are accomplished with sensitivity time control and a vertical antenna 
pattern designed for uniform illumination as a function of depression angle. 

Deviations (from the ideal) of these functions can cause radiometric errors. 

b. Sensitivity Time Control (STC) 

Reference T contains instructions on how to adjust the STC to give the desired 
signal. This description includes wave shapes and is considered the best data 
available. The STC so adjusted requires no correction. The STC is turned off in 
modes used above 30, 000-ft altitude and does not apply to the flight of 30 August 1973. 


a USAF T.O. 12P3-2APQ102-2-4. Radar Mapping Set, AN/APQ-102 and AN/APQ-102A 
(Frequency Converter-Transmitter CV 1678/APQ102), Chapter 11. 
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e. Antenna Uluminat ion 


It can be shown that if the vertical antenna pattern has a gain 

2 l/2 

G = K esc 9 cos 9 , (29) 

then the terrain would be uniformly illuminated as a function of depression angle 9. 

Such a pattern can be synthesized over a limited angle. In the AN/APQ-102A, the 
pattern is normalized at 18 deg. Figure 4 shows the theoretical vertical antenna 
pattern of the AN/APQ-102A, together with the tolerances in gain. The standard 
deviation of the one-way gain from uniform illumination is less than 0.5 dB. The 

antenna pattern of the arrays to be used can be measured and the deviation between 

2 l/2 

measured values and the ideal csc“ 9 cos ' “ 9 determined. It was anticipated that 
the measured antenna pattern could be used to make radiometric corrections to the 
imagery of Wallops Island made on 30 August 1973. However, it has been deter- 
mined that USAF records do not make this possible. Therefore, no corrections for 
the antenna pattern were made. The image distortion program discussed in 
Section III has such provisions. 


1 


V 


V. 

l 
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SECTION III - DIGITAL PROCESSING OE IMAGERY 


1. DIGITAL AZIMUTH PROCESSING 

The processing described in this section will accommodate AN/APQ-I02 radar data 
which has been range compressed, recorded optically, scanned, and digitized for 
processing. The azimuth compression will be performed by a high-speed digital com- 
puter. The data flow diagram is given in Figure 5. The processing will be performed 
to obtain 30-ft-resolution imagery with the option of either one or two azimuth looks. 

If more rapid processing of the data is desired, the azimuth resolution may be degraded. 

The data will be processed for an azimuth offset that is nominally PRF/4. In actuality, 
the azimuth offset frequency is not exactly known, and the data processing must take 
this into account. To maintain low sidelobe levels, an oversampling factor of at least 
four will always be maintained. The sampling rate of the input data will be not be 
reduced until azimuth compression is being performed. 


The implications of the foregoing may be better understood by examining the fundamental 

formulas for azimuth compression. The synthetic aperture length which must be flown 

to attain a desired 3-dB azimuth resolution, \V. , is 

A 


0. 88AR 


L 


SYN 


2VV 


(30) 


A 


where R is the slant range to the target measured on a line perpendicular to the flight- 
path, and X is the wavelength of the transmitted signal. 


However, when the phase history of a point target which has been collected over the 
required is compressed (i.e. , processed in a matched filter), the resultant side- 

lobes of the sin (x)/x compressed waveform have a -13.6-dB peak and decay slowly. 

PAGE blank not hlmeh 
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Figure 5 - Digital Processing of SAR Data-Flow Diagram (Sheet 1 of 2) 


- 24 - 













SECTION III 


GERA -2 089 


0 


INPUT DATA 
DECK 


GROUND RANGE 
RADAR DISTORTION 
ANALYSIS PROGRAM 


PRINTED 


TEST TAPE 

OUTPUT 


GENERATOR 





! 



i 


IMAGE DISTORTION 

f CARD 




INPUT 


CORRECTION 


PROGRAM 







ORIGINAL PAGE IS 

or poor quality 


Figure 5 - Digital Processing of SAR Data-Flow Diagram (Sheet 2 of 2) 
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These sidelobes present a problem in that the sidelobes of a large target may have 
larger amplitudes than the mainlobe of a target, and hence mask it. To prevent this 
from occurring, a weighting function is applied to the return phase history. Basically, 
the weighting function reduces the sidelobes by the application of a symmetrical taper 
across the azimuth phase history of a target. This symmetrical attenuation, however, 
causes a broadening of the mainlobe of the target when azimuth compression is per- 
formed. To maintain the desired resolution, yet achieve a low sidelobe level, excess 
azimuth bandwidth is required. 


Typical of the weighting functions which may be applied are Taylor aperture functions. 
A Taylor aperture function which suppresses the peak sidelobe to -30 dB will broaden 
a point target's mainlobe by a factor of 1.42. When this weighting function is used. 


= 1.25XR 
L SYN 2W 

A 


(31) 


Processing is being performed for a finer resolution than is specified, with the know- 
ledge that the weighting function is being utilized and will broaden the mainlobe to that 
desired. 


If the antenna's real azimuth beamwidth, ( 3 , (here considered to be 3-dB beamwidth) is 
capable of illuminating more azimuth extent than is required for a desired resolution, 
i_.e. , if 


R (3 > I 
s 


J SYN 


(32) 


then the excess illuminated area may be used to form more than one synthetic aperture. 
For the AN/APQ-102A system, the 2-way 3-dB beamwidth is approximately 1 degree, 
which would theoretically imply that 8 synthetic apertures for 30-foot resolution exist 
between the 3-dB points of the antenna beam: 


R /? R (1 deg) / KXR 
s s / s 

L SYN = 57 - 3 / 2W A 


8.37 


(33) 
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However, because the limitations of the motion compensation INS, azimuth recorder 
bandwidth, etc., blurring of the image will possibly occur if more than two looks are 
combined. Thus, for this problem, 

L_. rvr = 2.083 X 10 3 X R x number of looks , (34) 

SYN s 

where L is the synthetic aperture length in feet, R is the slant range in feet, and 

ol JN 

the number of looks is either one or two. 

The AN/APQ-102 radar has a PRF of 1. 1 V, where V is the aircraft velocity in feet per 
second. Therefore, a sample of the terrain is collected once per 0.9091 ft of aircraft 
travel. This sampling is greatly in excess of that necessary for 30-ft resolution, which 
is 

30 

minimal sample spacing = — — — — ft . (35) 

2 X 1. 25 X number of looks ' 

(The factor of 1. 25 in the denominator accounts for the excess bandwidth required to 
preserve resolution when using the weighting function.) It is necessary to have this high 
PRF to keep the spectrum which lies within the antenna's mainbeam unambiguous. 

The unambiguous bandwidth of the sampled spectrum lies from zero frequency to l/2 
PRF, or 0.55 cycle/ft. The clutterlock, however, keeps the antenna, and hence the 
doppler spectrum within its mainlobe, centered on zero doppler. Therefore, it is neces- 
sary to translate the return mainlobe spectrum up in frequency so that a frequency of 
-0.275 cycle/ft will lie at zero frequency, and a frequency of ^0.275 cycle/ft will lie at 
0.55 cycle/ft. This is accomplished by mixing the return with an azimuth offset fre- 
quency of l/4 PRF. The translated return spectrum is illustrated in Figure 5. Energy 
at frequencies above ±0.275 cycle/ft will fold back into the spectrum of interest. 

However, as is illustrated by the dashed lines in Figure 6, this energy is heavily 
attenuated by the rolloff of the antenna's mainlobe. The peak sidelobes of the AN/APQ- 
102A 2-way antenna pattern are more than 26 dB down, and thus the energy in them will 
contribute little to the processing noise. 
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Figure 6 - Antenna Azimuth Frequency Response 

The spatial bandwidth required to process 1 look for 30-ft, 3-dB resolution with 30-dB 
peak sidelobes is (1.42)(0. 88)/(30) = 0.04617 cycle/ft. The factor of 1.42 is the 
excess bandwidth ratio required for the sidelobe control. Thus, the bandwidth of the 
spectrum which will enter the processor is in excess of that necessary to process one 
look by the ratio of (0. 55)/(0. 04167) = 13.2. As will be detailed later, the input band- 
width can be reduced by filtering to reduce the excess bandwidth and improve the signal- 
to-noise ratio. (Note that only 8.37/13.2 of the unambiguously sampled data lies within 
the 3-dB antenna beamwidth, as was shown in Equation (33)). 

The radar’s azimuth bandwidth and sampling rates having been examined, the effects 
of scanning the optically recorded data and digitizing samples will be considered. For 
each range element, an azimuth sample has been taken and recorded on film for each 
0.9091 ft of aircraft travel. When this data is scanned, however, a sample is taken from 
the film at the equivalent of once for each 0. 826 ft of aircraft travel, which introduces 
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an effective increase in the azimuth sampling of a factor of 1. 10. It must be understood 
that no increase in the information bandwidth has occurred, that having been restricted 
by the original PRF. However, there is a translation of all frequencies because of the 
resampling. All data must be treated as though the spatial bandwidth were (0. 55)(1. 1) = 

0.605 cycle/ft, even though no information lies in the portion of the unambiguous spec- 
trum resulting from the different input and output data rates. 

The digitized data will be treated as if the original PRF produced a factor of (13.2)(1. 1) = 
14.52 in excess bandwidth over that required for a single azimuth look. Each azimuth 
look will thus occupy (0.04167/0.605) X 100 = 6.9 percent of the unambiguous sampled 
bandwidth. 

The first operation in the digital processing of the data is to bandpass only those 
frequencies necessary for azimuth compression, and thereby improve the signal -to-noise 
ratio by reducing the noise bandwidth, 'piis is done in the azimuth prefilter, which has 
been designed (described below) to have a nominal center frequency of 0.4545 (i.e. , 
PRF/ 4 / 1 . 1) of the sampled bandwidth and have 0. 04167 and 0. 08334 cycle/ft spatial 
bandwidths for the one- and two-look cases, respectively. 

3 . 

The two azimuth prefilter functions are "window function" designs. To produce a 
window function with a desired center frequency and bandwidth, the following procedure 
is employed (the steps are illustrated in Figure 7): 

1. in the frequency domain, place an impulse at the desired fraction of 
the bandwidth for which the filter's center frequency is to lie, and at 
the corresponding negative frequency 

2. Take the inverse discrete Fourier transform (ID FT) of this spectrum. 

The result will be a sampled cosine in the time domain having an 
integer number of cycles over the time extent of the IDFT output 


Gold, B. and Rader, C.M.: Digital Processing of Signals . McGraw-Hill, 1969, pp. 217-231. 
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3. A weighting function, in this case, a 40-dB Taylor function, is point 
by point multiplied with the sampled cosine. This is the filter function 

4. A discrete Fourier transform (DFT) is taken of the product. This 
shows the bandpass of the filter. Note that if no window function were 
applied, the step 4 output would have the shape of sin (x)/x, which is 
the transform of a pulsed cosine. 

The selection of the number of points in the window function design depends essentially 
upon three factors: 

1. The accuracy to which the center frequency of the filter must be 
positioned 

2. The bandwidth of the filter 

3. The rolloff rate and minimum stop-band attenuation of the filter. 

For the two azimuth prefilter functions, the desired center frequency of both is 0.4545 = 
5/ll. Thus, as an impulse is needed at both positive and negative frequency, the length 
of the reference function should be a multiple of 2 X 11 = 22 points. 

The bandwidths desired are 0.04167 and 0.08334 cycle/ft spatial bandwidths. As the 
3-dB width of the sin (x)/x function mentioned in step 4 of the foregoing is 0.88 X 2/N of 
the bandwidth, where N is the number of points in the sampled cosine, and as the 40-dB 
Taylor weighting broadens this by a factor of 1.42, then for the 0. 04107 cycle/ft filter, 

0.88 X 2 X 1,42 tN = = 40 , (36) 

and for the 0. 08334 cycle/ft filter, 

0.88 X 2 X 1.42 tN = ^-N= 20 . (37) 
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Thus, for the filters, N = 44 and N = 22 provide excellent choices, with the impulses 
at (10, 34) for the first and at (5, 17) for the second. The taper over the amplitude of 
the filter's output spectrum will be used as part of the weighting to control the azimuth 
sidelobes (described in Equation (31)). 

The two filters are readily translatable to different center frequencies simply by changing 
the position of the impulses in step 1. The narrow bandwidth filter can be stepped in 
increments of 0.04545 of the sampled bandwidth and the wider bandwidth filter in steps 
of 0.9091 of the sampled bandwidth. Finer steps can be obtained by increasing the num- 
ber of points in the reference function and changing the window function accordingly. 

The azimuth prefilter is a nonrecursive, convolution filter. In nonrecursive filters, the 
output is not fed back to the input. Although filters with feedback (i_.e. , recursive filters) 
have shorter reference functions than do nonrecursive filters, they suffer in that they 
allow a noise buildup because of signal quantization and do not offer the truly linear phase 
characteristic which nonrecursive filters provide. Hence, nonrecursive designs are 
considered superior for this application. 

In the nonrecursive filter, N consecutive data points are multiplied by the N corresponding 
filter reference function points; the N products are summed, and the result is obtained. 
The oldest input data point is discarded, the remaining N - 1 data points are shifted one 
position, a new data point is entered, and the multiplication and summation process is 
repeated. Thus, for every data point entered, there is one data point output. 

After azimuth prefiltering has been performed, the data will be compressed to its ulti- 
mate azimuth resolution. The length of the synthetic aperture required to compress an 

_3 

azimuth return was shown in Equation (34) to be equal to 2.083 X 10 X R x number of 

s 

looks, where R is the slant range to the target in feet. A digitized data sample is taken 
s 

from the data film once for every 0. 826 ft of aircraft travel. The number of data points 
contained in the synthetic aperture length, N, is given by 
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N = 2.52 X 10 ** X li^ X number of looks , (38) 

where R is in feet, 
s 

The azimuth compression filter function with which the prefiltered data will be correlated 

will next be determined. Consider an isolated point target at a slant range, R . The 

s 

ground range to this target, R , is 

S 

R = (R + 32.8 X M) ft , (39) 

g v go 

where R is the ground range (in feet) to the near edge of the swath being mapped, 32. 8 
is the conversion factor from film scan to feet on the ground, and M is the number of the 
range cell in which the target lies (M equaling zero for the first range cell). The slant 
range and ground range are related by the equation 

I 2 2 1 1/2 

R s = | h gl ' (40) 

where h is the aircraft altitude. 


As the aircraft flies past the point target, the phase of the return, 0 , is equal to 


4ttR 

<t> - , (41) 


where R is the slant range to the target, and \ is the radar wavelength. The range R 
s 

may be expressed as 



(42) 


where R gQ is the slant range to the target when measured on a line perpendicular to the 
flightpath (Ue. , at the point closest to the aircraft), and X is the along-track displacement 
of the aircraft from the point at which R g is measured. As R gQ >> X, Equation (42) may 
be approximated as 
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R * R + 
s so 


2R 

so 


with a high degree of accuracy. Hence, 


i 4tt 

*■ = r 


R 


x 2 


so 2R 


so. 


2ttX 

= 6 + 

o AR 

so 


(43) 


(44) 


where 6 is a constant, 
o 

To Equation (44), the azimuth offset frequency (shown previously to be PRF/4/l. 1, where 
the factor of 1.1 results from the digitizing process) has been added. Thus, the azimuth 
phase history of the signal presented to the azimuth compression filter is 



(0. 82 6n) 


2 



+ 32. 8M 



1/2 


- 2tr 


(0.227n) + 0 


(45) 


where 

n = azimuth sample number, n = 0 occurring at X = 0 
0.227 = azimuth offset frequency in cycles per foot on the ground after digitizing 
0. 826 = distance between samples in feet. 

The azimuth compression reference function (ACF) which will compress the point target 
to the desired resolution is 

ACF = exp | j^| = cos 4> + j sin 4 . (46) 

The value of 4> o for the reference is set to zero as it is an arbitrary constant. The value 
of N has been determined in Equation (38). The value n in Equation (45) will be stepped 
from -n/2 to N/2 for generation of the reference function. The compression is performed 
by the complex convolution of the data and the reference function, although the data quad- 
rature component is always zero, and hence no multiplication is performed with this term. 
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The computed ACF will be weighted by a Taylor aperture function to reduce the azimuth 
sidelobes. (Recall that the broadening of the mainlobe of the compressed pulse has 
already been compensated for by the factor of 1.42 in the aperture length formula. 
Equation (31).) The Taylor aperture function has only real, positive coefficients. 

The product of the weighting function and the ACF will result in a function 

ACF weighted = A<n) exp M • < 47) 

where A (n) = 1 when n = 0, -N/2 « n e N/2. 

max 

The output of the azimuth compression convolution will be generated having 7. 5— ft 
spacings, or equivalently at one-ninth of the input data rate. This will reduce the data 
rate and hence the number of calculations by a factor of nine, yet retains a sufficient 
number of data samples to preserve the processed resolution after detection. 

Detection of the compressed data is accomplished by forming I 2 + Q 2 of the azimuth 
processed image; i.e. , by squaring the real and quadrature components of the 
data and then summing them. Detection produces information which contains only 
magnitude information, the magnitude being proportional to the power of the return over 
the aperture length from a point target. 

To obtain two looks in azimuth, the ACF will be twice the length as that used for one- 
look processing. The weighting function is applied in the same manner; however, twice 
as many sample functions are taken over the Taylor aperture function. The two looks 
are formed after azimuth compression and detection by passing the data through a post- 
detection filter. The postdetection filter's impulse response is equivalent to 30-ft 
resolution. This filter is formed by summing four consecutive azimuth compression 
filter outputs and dividing by four; i_.e. , 

4 

E 0ut - 0.25 £ E(l> . (48) 

1=1 

For every output of the azimuth compression filter, there will be one output of the post- 
detection filter. 
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For the azimuth compression, it has been shown that the phase of the ACE varies with 
range as shown in Equation (45), and that a new ACE can be computed for each range 
increment. Experience has shown, however, that satisfactory results may be achieved 
even if the phase difference between the signal and the ACF varies by as much as 
±22-1/2 degrees. Hence, 


<t> = 


4-ttR 


, 2ttX 

0 + XR 


Therefore, 


dci 2 -itx tr 

dR = i dR o = 4 
0 XR 


, KXR 

1 _ 0 

X " 2 SYN ' 4W . 


„ KXR 2, 

h. 2. AR _ _ 1 

2 4W AR " 2 AR " 4 

XR A 8W" 


2W 


lARl = 


11.520 


ft 


K^X (number of looks)' 1 ' 

For one look, AR = 57G0 ft; and for two looks, AR - 1 1440 ft. 


(49) 


2. REMOVAL OF IMAGE SKEW 

Because of such factors as aircraft across-track motion, antenna pointing errors, errors 
in scanning the data from the film, etc., the output data may be at a skew angle. The 
skew is corrected by the image distortion correction program. 
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The skew removal is accomplished by a "zero data" addition procedure. From the 
geometry shown in Figure 8, where p is the skew angle, it is seen that the data must be 
rotated to orient the data's range vector with the vector perpendicular to the flightpath. 
To accomplish this, data points with magnitudes of zero are inserted at near and far 
range to form a rectangular data block, as illustrated in Figure 9. (The all -zero 
columns at near and far range are for computational convenience.) 

The number of zeros added to each range return is given by the equation 

NZ = NRG sin (|$l) + 2 (50) 

where 

NZ = the number of zeros added for skew correction 
NRG = the number of ground range sweeps 
|$! = the absolute value of the skew angle. 

The skew correction is calculated by 

Y = AX ♦- B , (51) 

where 

Y = the corrected ground range 

X - the along-track distance relative to the patch being deskewed 
NZ 

A = ^* SIGN <*> 

_ 1 if SIGN (4>) is positive 

NZ if SIGN ($) is negative. 

Finally, it is noted that the number of zeros added to each range return must be an 
integer. If the number of zeros added to the ith range return is IZ., and the number of 
zeros added to the same return at far range is JZ., then 
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,z i ' l' A >«" 1 "luoilNUEn 

JZ. = NZ - IZ. 

1 1 

The placement of the zeros is illustrated in Figure 10. 

3. SCALE FACTOR CORRECTION 

Upon completion of the azimuth compression, the output data sample points may be 
spaced differently in range and azimuth. The image will consequently appear distorted 
because of the differing range and azimuth resolution. To compensate for this, a scale 
factor correction may be necessary. 

Scale factor correction is accomplished by linear interpolation on the azimuth data. For 
example, assume that an azimuth sample was calculated every 12 feet, and that 30 feet 
was desired between samples in both dimensions. Then, to achieve azimuth samples 
spaced by 30 feet. 



(52) 

(53) 
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f 



T 



-(Y + Y > 

4' 3 4 ; 


f 



etc. 


(54) 


It is observed that data points Y and Y are not utilized in the foregoing calculation. 

A D 

Therefore, an increase in the processing rate is possible, because these points need 
not be calculated. 
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SECTION IV - CONCLUSIONS 

The results of the geometric distortion analysis indicate that the distortions in AN/APQ-102A 
imagery are primarily the result of navigation system errors that are external to the radar 
system itself. These distortions can be rather high in magnitude (e.g. , one percent), but 
have a low spatial frequency. As such, it is a relatively simple task to measure and remove 
the geometric distortion. Measurement is accomplished by comparing image distances 
(obtained from a map) between known ground features with good distances. Using this tech- 
nique, the residual distortions were under 100 meters. Computer programs to measure and 
correct these distortions were delivered as part of the contract effort. 

The major program effort consisted of generating a computer program to digitally process 
AN/APQ-102A phase history data. This program was checked out and validated at the 
customer's facility — thus providing a capability of processing subsequent AN/APQ-102A 
data without contractor support. 
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C P R 0 G • 9 - 7 5 POINT TARGET GFNERATOR 

C 

DIMENSION IDATA ( 8800 )#MTR1(27) 

C 

S C N * • ? 

BIS. .316 

ALGF» .A3A?94481?03?? 

PI *3. 1415926535897 
Al AN. . 10?? 
t}V.»5*PI/180. 
fpal«a. *pi/alam 
REWIND 2 

c 

200 CONTINUE 

RE AC 98 » MR 1 
PRINT 1C1 
PRIM 102 

PRINT 99/MTR1 . 

C OUTPUT RECORD 1 ON TAPE 

C 

CALL TAPE R1 (MR1 ) 

C 

C 

C CAZ DISTANCE per AZIMUTH SAMPLE 

C DRG CISTANCE FFR RANGE SAMPLE 

C A Z0FS azimuth OFFSET 

c 

C HF ALTITUDE IN FECT 

C RO RANGE IN FEET 

C PHIO FHASE ANGLE 

C 

READ 100/DAZ,DRG/ AZOFS 
PRINT 103/DAZ/DRG. AZ0FS 
READ 1CC*HF,P«,PHJP 
PRINT lO*/HF,Ke,PHI0 

201 CONTINUE 
C 

C M RG ELEMENT 

C N AZ ELEMENT 


ORIGINAL PAGE IS 
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C NRe NUMBER OF RANGE’S 

C 

RFAD 105# M#N# N'RP 
I F ( NRB ) 198# 198# 199 

198 CONTINUF 
CALL EXIT 

199 CONTINUE 

PRINT 106#M#K-#NRB 

I X » 1 

XM»M 

XN«N 

C RG LOOP 

C 

D« P02 I y « 1 # N'K'B 
CO 2C3 1 ■ 1 # 880C 
IDaTa ( 1 ) «C 
203 CONTINUE 

c calculate az index limits 

c 

fiS»SCKT(MF*MF*(Re+XH # DRG)«*2) 

NA2-SIN(PW)*RS/CAZ 

NA Z 1 »N-NAZ 
\AZ2«N+NAZ 
C Ll^lT INDICES 

c 

NAZ1*NAX(1#NAZ1 ) 

N A Z 2 ■ M I n ( N a Z 2 # 17600) 

PRINT 107#NAZ1#N AZ2 

C ENC OF BOOK KEEPING# D0 THE CALCULATIONS 

C 

wAZ«NAZl 
A N ■ - N A Z 

AA».5»FPAL»DA?*DAZ/RS 
BB»-2.*PI»AZ0FS 
PRINT iC9# AA # B13 
IY«C 

IF (NAZI -8800) 208# 209# ?09 
C OUTPUT FIRST HALF 0F AZ SWEEP 

C 
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2 09 CONTINUE 

CALL TAPEL2 ( IDATA, IX/ I V ) 

I Y»?8C0 
JA2«\AZ1-8800 
208 CONTINUE 
C AZ LOOP 

C 

DO 20A JAZ.VAZWNAZ2 
PHI * ( AA*AN>PB )*AN + PMIP 
C TAKF CSS , ADD BUS AND SCALE 

C 

T=SCN«COS(PHI ) ♦? I S 

c convert transmission to density 

c 

D=ALGP*ALPG( l./T) 

C NOW scale t« 8 bits 

c 

1D*D* l28»+»5 
I Data ( jaZ ) * I D 
c 

GO TO 197 

PRINT 1c*#AN#PHI#T#D# ID 
197 CONTINUE 
A N « A \ + 1 • 

JAZ-JAZ^I 

C 

IMwA7-88C0)20A#204#205 

205 CONTINUE 

C OUTPUT MALE OF AZ SWEEP 
C 

CALL TAPEW2( IDATA# IX# IY) 

DO 2q 6 I ■ 1 #88oo 
I DATA ( I ) -0 

206 CONTINUE 
IY-S8C0 
JAZ ■ 1 

c 

20A CONTINUE 
210 CONTINUE 


is 
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C auTPUT SECOND HALF OF A Z SWEEP 
C 

CALL TAPEW2( IDATA# IX/ IV) 

IF ( IY >?lc/2l2/2U 

212 CONTINUE 
I Y * 88GO 

CP 213 1*1/8800 

I P A T A ( I ) .0 

213 CONTINUE 
GP TP 210 

£11 CONTINUE 
C 

I>*lx*l 

x^-x^+l . 

202 Continue 
C 

GP TP 201 

98 fph*aT(8CAI) 

99 FPRVAT ( 6X80A 1 ) 

IOC FOPNAT { 8E 10 • A ) 

101 FPR^AT (6X1 AWPK0GRAM 9-75/ 2OX3OHG0ODYE AR AEROSPACE CORPORATION#/) 

102 FPRNAT(6X22HPeiNT TARGET GENERATOR#/) 

103 FORMAT ( 6X5HP A £ *El 2 • 5/ 5X5HDRG *E 1 2 • 5, 5X5HAZ0FS/ El2 *5# / ) 

10A FyRNAT (6 x5hmF -F 1 2 • 1 / 5x5hR0 »F12. 1#5X5hPhI0»Fi2.3/ / ) 

105 FPRNATU6I5) 

106 FPRNAT(6X,3HN .I7/5X/3HN « 1 7 / 5X3HNRB# 1 7, / ) 

1 0? FPRPATI6X21MAZIPUTH SAMPLE L I MI TS/ 2 1 1 0# / > 

108 FPRNAT(AF12.A, 15) 

109 FORMAT (6XAHAA .£12.6/ 6XAHBB *E12»6#/) 

END 
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SUt3R8UT I NE TAPEPl(HTRl) 

DIMENSION MTR1 (27) 

THIS ROUTINE IS FOR WRITING RECORD 1 W I S S MAG TAPE 
FORMAT (27 CHARACTERS) 

27 CHARACTERS 0F WHICH 20 OR 22 ARE NEEDED 

"RITE TAPE 2/ MTR1 

RETURN 

END 


SUBROUTINE TAPEW2( ITATA/ IX, IY) 

C THIS R9UTINE IS F OR OUTPUT TO MAGNETIC 

C TAPr FOR THf SECOND AND SUBSEQUENT RECORDS 

DIMENSION IDATA (8800) 

WRITE TAPE 2, IDATA 

return 

100 FORMAT( 16, IX, 4013) 

RETURN 

END 
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AZIMUTH PROCESSING PROGRAM 
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Card Input to the Azimuth Processing Program 

Card I 80 column description card. 

Card 2 N number of data points/record. 

(2 values) (Azimuth samples per range bin.) 

NRG number of range bins to be processed. 


Card 3 SCL * density represented by a value of 255. 

(3 values) TMEAN » mean transmission for bipolar calculations. 

SCL2 ■ second scale to convert to integer. 

These are for Initialization of the Density to 
Transmission to Bipolar Routine. 

Typical Data: SCL = 2 or 3, TMEAN = . 38 , SCL2 = 61* or 128. 


Card A NL00K ® number of looks, l or 2. 

(3 values) K = K in COS (2P I K l/N], (see azimuth prefilter formula). 

NB - number of bits in the quantized output. 

Card 5 PHlO » phase offset in degrees. (This is an arbitrary input.) 

(1 value) 

Cards A and 5 are for the azimuth prefilter reference function 
generator routine. 


Card 6 DAZ = distance per azimuth sample in feet. 

(3 values) DRG = distance per range sample in feet. 

AZOFS = azimuth offset frequency in cycles/ft on ground. 


Card 7 HF = altitude in feet. 

(3 values) R0 = range in feet. 

PHIO = phase angle in radians. 


Cards 6 and 7 are for initialization of the compression 
reference generator routine. 
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C MAIn LInE PROCESSING or imagery 

c 

D I MENS I ON IR<44># JD { 45 ) 

DIMENSION NM(20>* RR ( 44 ) 

C 

dimension in<i7f00) 

DIMENSION KO (931) 

DIMENSION I CR ( 930 ) / ICQO30) 

C 

EQUIVALENCE ( ID(1>#JD<44> ) 
EDUIVALENCE (KC(931 )# I D ( 1 ) ) 

N = 1760C 
jT *1 
I T a? 

C TAPF ? IN, TAPE 1 OUT 

2 00 CONTINUE 
REWIND IT 
REWIND jT 

c nn; m >*q col discretion 

READ 99# NN 
PRINT 1 cl 
PRINT 9«#NN 
READ 100* ni#T'RG 
PRINT 102#N#NR'G 
NT »N 
< = C 

CALL DNTRBP <1D#K) 

CALL REPGEN (IR#NPR) 

C 

NLOOK ■ 1 

IF(\P«-?2) ?20# 220# 221 

220 CONTINUE 
NLO0K • 2 

221 CONTINUE 

M • -NLPOK 

CALL CMPREP ( I CR# I CQ# M ) 

LOOP ON RG 
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pax»o 

cp 201 KRG*1/NRG 

N » KT 

CALL T'APE IN ( I D> NT ) 

K ■ N 

CALL DNTRRP ( ID/<) 

CALL PRINTD { I 0 > N ) 

C 

NP ■ \ 

CALL AZPRF ( ID# TR*N#NPR»JD/NP) 

SHIFT Data by 9 BITS 

DB 207 I a \ , n f ** 

J0( ! ) « wD(I >/512 
207 CONTINUE 

call printd (JD>N p ) 

CALL PL(?TD( JD#NP) 

C 

P a *PG 

Call CMORrF ( ICR/ ICO/*) 

\R « P 

CALL PR INTO ( I CR/ M > f 

CALL FRINTO (ICG/^) 

C 

NS ■ 9 

CALL AZCSMP ( jd, ICR# ICQ/NP/NR/NS/KD) 

\ a NP/NS 

GP TP (224,222), NLP0K 
222 CONTINUE 
C 

C TWO LP3K CALCULATI0NS 

C 

\ P ■ N • 3 

OP 223 I ■ 1# NP 

«C(I> ■ KC( !)*KD( 1+1 )+KD( I *2) ♦K0< 1+3) 

* 223 continue 

In ■ NP 

224 C6NTINUE 

li 
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CALL PRTMD (KD<N) 

D p cC6. I«1#V 
IF ( r Ax - kD ( I ) >205*206/206 
c05 CONTINUE 
MAX«KD(I) 

PRINT 10^/NRG# I/MAX 
cOt> continue 

\P*\+1 

?03 I «NP, NT 
<F < I ) s 1 
203 CPNT1NJE 

CA(_l TapEPUT (<D/NT ) 
cOl CONTINUE 

CALL PLOTD (KC,N) 

Gfi TO 2oC 

CALL PL0TC ( ID/K) 

aRITE TAFE JT, C JO { I ) # I • 1 # < ) 

99 FORM AT (6X20AU,/) 

99 FORVAT (20AA) 

100 FfiRyAT (PF1C*4) 

101 F&RMAT ( 36X35H I ^ AGE PROCESSING G»A#C« PROG 9-75///) 

102 FORMAT (f,Xl9MNUHRER OF SAMPLES ■ * 1 7# 5X22HNUMBER OF RANGE BINS »I7 
1 #/) 

10 4 * FPRPAT <U0*2xlCP) 

END 
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SUBROUTINE DNTRPP ( I D # M ) 

DENSITY TO TRANSMISSION TP BI PPLaR 

DIMENSION' ID ( I7f CO ) 

DIMENSION JD( 256 > 

I E ( M ) 2oC/ 200* 210 

200 CONTINUE 

scl » density represented by a value of 255 

TM£AN « m EAN TRANSMISSION FOR BIPOLAR CALCULATIONS 

SCL2 « SECOND SCALE T8 CONVERT TO INTEGER 

TYFICAL Data: SCL «2 or 3/ TMEAN-.38# SCL2*64 OR 128 

READ 100* SCL# TMEAN# SCL2 
PRINT lol 

PRINT lu2# SCL# TME AN# SCL2 
PRINT io3. 

SCL • SCL/256 • 

DP ?C 1 I ■ 1# 2P6 
< » 1-1 
X « < 

X » X * SCL 
X » 1 C • * * ( -X ) 

X « X-TmEAN 
JD( I ) - X • SCL? 

PRINT ICA* «» JO <1 ) 

201 CONTINUE 

RETURN 

210 CONTINUE 

\ a M 

BIAS « 0 # 

DO 211 I - 1# N 
< » I D ( I ) ♦ 1 
I E ( K ) 212# ? 1 2 # ?13 
213 CONTINUE 

IE ( K-255 ) 214# 214# 212 
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212 CONTINUE 

PRINT 105/ I0( I ) 

214 CONTINUE 

ID( I ) ■ JD(K ) 

X * ID< I ) 

BIAS » III AS*X 
211 CONTINUE 
X « N 

BIAS » FlAS/V 
PRINT 106/ c I AS 
RETURN 

100 E*R!*AT (8E1C.4) 

101 FORMAT (/////, 6X40H0ENSITY T8 TRANS^ISSI 8N T8 Bl PBLAR DATA,/) 

102 FORMAT (6y,5MSCL « / Flo* 4# 5X5HTKE AN# Flo* 4#5X5HSCLH»* FlO* 4# / ) 

103 FPRyAT (25X7HDATA IN, 10X6MRESULT# / ) 

104 FORMAT (25XI5/1CXI5) 

105 F0RPAT (6X33H**. INPUT DATA 8UT 8F RANGE ***#15) 

106 FORMAT (6X31HAVERAGE RlAS F8R THIS RANGE WAS/FlO*4) 

END 
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subroutine: bergen (iR/N) 

AZIMUTH PREFILTER REFERENCE FUNCTION GENERATOR 
TWO DATA CARDS REQUIRED 

DIMENSION IR(AA) 

DIMENSION RR(AA) 


nlopk ■ number *f looks 

K ■ K IN COS ( 2P T K I/N) 

\r ■ NUMBER of BITS IN THE CUANTIZED OUTPUT 

PHIO > PHASE OFFSET 


PI - 3* 1^15^26535897 
READ 103/ NLOPK/ K/ NB 
°p I NT 1 0 A / N'LppW/ K/ NB 
READ 100* PH 1 0 
PRINT 105/ PHIO 
SC « ( 2 • **NE ) • 1 • 

N » 2? 

Ge TP (201/ 20? ) / NLO0K 
201 CONTINUE 

N a A 4 

20? CONTINUE 
XN • \ 

PRINT \ c e 
X^IN 3 C» 

XMAy a \ -1 
A = AC. 

NB a 6 
X » C • 

CALL TAyLPR ( XM ! N / X^AX/ A/ NB/ X/ WTN/ 0 ) 
DP 203 I « 1/ N 
A a < I • 1 ) *< 

A a 2 • *P I * A/XN + FH I 0 
X a 1-1 

Call TAYLOR (XMIN/XMAX/A/NB/X/WTN/1 ) 
RR( I ) a WTN*C0S(A) 


sssfigsf 
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203 CONTINUE 

PUT amplitude wfirhting here 

2C A I « 1 » N 
F5 « «5 

I F ( PR ( r ) ) 205/ 205/ 206 

205 CONTINUE 
P5 ■ “ • ^ 

206 CONTINUE 

IP t I ) » Rp(! )«Sr*P5 
print 107/1, IR(I) 

20* CONTINUE 

return 

100 FORMAT ( P E 1 C • * ) 

103 FORMAT (1615) 

10* EpRRAT (6X7HNL6PK -/I3/5X3HK ./I3/5X*HNB ./I3//J 

105 FORMAT (6X7HPWI? »/Fl0»*//> 

106 F5R^AT(//,6X9HRFF.ERENCE#10X'1HI,10X1HR,/) 

107 FflRKAT (-15X2-1 1 1 )' 

END 
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SUBROUTINE AZPRT ( 1 0* I R # N, NPR, JD# NP ) 
C A7IVUTH PREPILTFR WALLBPS DATA 

DIMENSION I p { 4 4 ) / JO ( 45 ) 

C I I'.pNG I BN I D ( 17600 ) 

C 

. NP»N-NPR 

CP 2C? !«1#K’P 
IS«C 

oe ?C3 J»1/A'PH 
K « I+J-l 
IS«ID(K)*IR(J)*JS 
c03 CONTINUE 
JD( T >*IS 
c02 CPNTINUP 
RETURN 
END 
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C 

C 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBR6UTIKF CMPRFF (IR,IQ/M) 

COMPRESSION REFERENCE GENERAT9R 

DIMENSION IP(«>3C ), IQ (930 

M ■ -1 f)R -? FOR INITIALIZATION ABS<M) ■ NL0QKS 

m returned as number of points in the ref 

IF(m) 2oC# 21C * 210 
200 CONTINUE 

PI » 3«1415P26535SP7 
Al AM « . 102? 

'fi A « 3 0 • 

SNF « .88 
TWF ■ 1.4? 

XK * SNF«TWF 
XLOOK » IABS(M) 

PPAL * A.»PT/ALAM 


Daz « DISTANCE per AZIMUTH sample 

DPG ■ DISTANCF p ER RANGE SAMPLE 

AZOPS « AZIMUTH OFFSET 

HF . ALTITUDE IN FEET 

R? « RANGE IN FEET 

pm io « phase angle 


ND3C = D5 DPWN FOR TAYLOR WEIGHTING (40) 

NBaR ■ \BAR FOR TAYL0R WEIGHTING (6) 

NF-ITS ■ nUMSEk of BITS IN REF FUNCTIONS 6 * 7*8 

READ 100# DA Z, CRG, AZ8FS 
PRIM 103* DAZ# DRG, AZ8FS 
RE AC 100' HF, Re, PHI0 
PRINT 104, HF, Re, PH I e 
ND3C » 40 
NPAR ■ 6 
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Nr- ITS - P 

P P I M 105# N'DBD# NBAR# NBITS 
DPI) « \DPD 

SC « (2. **\'P ITS )-l . r 

PM * •5*Xk*ALAl^/(WA*DAZ ) 

P Nl o RN*XL0PK 
201 CONTINUE. 

RETURN 

2 1 C CONTINUE 
XN » y 

’S « SCWT ( HP #MF+ (R0*-XH*DRG >* *2 ) 

N A Z * RS«PN/2. + .5 
A f. » - N A Z 
N > NAZ+NAZ+l 
V « N 

X N I N 3 - N A 2 
XMAX = NAZ 

AA « •5«P p AL*nA7*DAZ/RS 
BP * -2 • *P I * A79FS 
P p J\T 106# AA# PB 

CALL TAvLPR (XNlN i XMAX#DBO#NBAR,AN, A mP#0) r 

CO 211 I « 1# N 

CALL TAYLPR (XPIN # X s 1AX#r)BD#NBAR#AN# AMP #1 ) 

PHJ « ( AA* AN+BB ) *AN + PHI9 
A \ * AN+ 1 • 

R « A^' p *SC*CeS(PHI ) 

0 » A y P*SC*SIN(PHI ) 

I f ( I) ■ R 

1 r. ( I ) ■ c 

211 CONTINUE 
RE-TLPN 

100 FPRKAT (firiO.H) 

101 FORMAT ( I10#AR1?.6#2I10) 

102 FPRN A T ( 1 6 I 5 ) 

103 FORMAT (6X5MDAZ » # E 1 2 • 5 # 5X5M0RG >#E 12 * 5# 5X5HAZ0FS# El2 • 5# / ) 

100 FORMAT (fcX5MHF • # F 1 2 . 1 # 5X5HR0 ■# F 12 . 1 # 5X5HPH If) ■# Fl2 . 3# / ) 

105 F9R y AT (6X6HNDB0 «, 1 5, 5X6HNBAR •# 1 5, 5X6HNBI TS 1 5, / ) 

106 FORMAT (6X4HAA ■# E 1 2 • 5/ 5X AHBB »#El2#5#/> 

END [' 
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SUBROUTINE A Z-CPMP (ID> IR# IQ/N,Np,NS* jd) 
c AZIMUTH COMPRESSION F I LTER 

c 

DIMENSION in ( 17600 ) 

DIMENSION JP ( 45 ) 

DIMENSION I R ( 93C > » 10(930) 

N « NUMBER OF POINTS IN THE DATA 
NR « NUMBER OF POINTS IN THE REFERENCE 

NS • NUMBER OF POINTS FROM OUTPUT SAMPLE TO NEXT OUTPUT SAMPLE 

ID IS DATA 

UP IS OUTPUT 

IR IS Real ref channel 

IG IS CUAD REF CHANNEL 

N * N-NR 

< a 1 

DO 20C I ■ 1# N, NS 
I SR ■ 0 
ISC » 0 

DO 201 J • 1 / NP 
L * J4 I • 1 

wR « I R ( J ) * ID ( L ) /256 
JC « I 2 ( U ) * I C ( L ) /256 
1SR ■ I SR* JR 
I S3 » ISOJO 
201 C^NTINJE 
<K « 128 
<K « 64 
1SG » ISC/KK 
I SR • ISR/KK 
IS « ISR*ISP*1S0*ISD 
UD(<)»IS 

< • K ♦ 1 

200 CONTINUE 
RETURN 
END 
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SUBROUTINE TAVLPR (XH! n ,xKAX,A,\B,X/AnS, J eB) 

THIS SUn^RUTINE CALCULATES THE TAVL0R APERTURE 

XH I K X W AX-X^IV r RANGF OF APERTURE 

XMAX 

A S DB DP WN PF FIRST SIDEL9BE IN T AYL9R ANTEnNa PATTERN 

NR = N-PAR (TAVLRR'S CONSTANT) 

X • value at WHICH 0NE POINT 6F APERTURE IS WANTED 
Ans ■ THE RETURNED ANSWER 

JOB « 0 T^IS IS A NEW XMlN#XHAX#A AND NB DATA SET 

* i usp previous xmjn;xmax#a and nb Data calculations 

DIKFnSIpN F l 2 C ) j C(10)# FF (2 ) 

EQUIVALENCE (FF(2)/F(1) ) 

IF ( W 9B ) 11# 'll# 12 
11 CONTINUE 

D a XHAX-X^IN 
PI • 3 « i ** 1 59 265 35897 
TP I « 2 • «F I 
^ ( 0 > ■ 

0? 10 I ■ 1# ?o 

A F G * 1 

F ( I ) = F ( 1-1 )*ARG 
10 CONTINUE 

a a * 10. ** ( aBS( a)/20. ) 

C a LL aR C®SH (AA,ETA) 

AA « ETA/PI 
A2*AA*AA 

ENB a NB 
NON ■ NP-1 

SIG » ENB/SQRt ( a2+(E n R-'5 )*«2) 

DO 20 N • 1* NBn 

FN ■ N 

EN2 ■ EN*FN 

S ■ 1 • 

DR 30 M • 1# NrH 

El'' - r 

ZM . SIG*SQPT( A2+(EH-.5)**2) 

S « S* ( 1 »-EN?/( ZM*Z M 5 ) 
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30 CONTINUE 

C ( N ) = S*E(rOM) # F(NBM)/<F{NBM*N)«r<NBM-N) > 

20 Continue 
12 continue 

F . (X-XNIN)/C*TPI-PI 
ANS ■ C» 

CMC N ■ 1, MBK' 

£N ■ N 

AN'S ■ A\S + C(N)«C0S(PN*P) 

40 CONTI Nut' 

AN'S ■ 1«4?««ANS 
«ETu«N 

PRINT ICC* A > NP 
PRINT 1C3* AA, SIG 
PRINT 101* V, C { N ) 

PPINT 102 

100 FPRNAT (5X17HSURR0UTINE TAYISR, //5X17HFIRST SIDEL0BE • *F 6»1# 
1 3H DE* l cX8N*N'-BAR • * 12///1OX12HC0EFFICIENTS) 

101 FORMAT (5X I5,5XF1S.8) 

102 FORMAT < 1 H 1 ) 

103 F0RNAT ( lByPHAA* 17 X 3HSIG*/2(5 xE15.8)#/) 

END 
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SUBR'UJT INE ARCPSH (U#A) 

UP = U.U 

!Faa-io ii/ ii/ l? 

12 CONTINUE 

A « AL3G(U + SCKT(U?-1« ) ) 

RETURN 
11 CONTINUE 

A , 0. 

PRINT 101 
RETURN 

101 FORMAT ( 1H0/5X56W*** TAKING ARCPSH fJF A NUMBER LESS THAN 1 I MP0SS I 
ISLE *•*»/) 

END 
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1 

2 
3 
* 
9 
6 

7 

8 
9 

10 

11 

12 

13 

1* 

19 

16 

17 

18 

19 

20 
21 
22 
23 
2 * 

25 

26 

27 

28 

29 

30 

31 

32 

33 
3 *. 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 
90 

51 

52 

53 


19135 APR 03# ' 75 
DOjBlE PRECISION £x<512)#EY(512> 
di-ension ibufiijqoi 

OOjBLE PRECISION XI 512 )#r I 51 2 ) #P I 512 ) > .. I SI 2 »# A# B# C# E I 512 > # D# F# 0# xl 
1(512>#"#aL#A1#a2#a3#a4# V#V1#A3I20)#XE(512)#YE(512)#C1#Y1( 

25121 

DI-EnSIOn aa< 16 l « **1 ( 20 1 # **21 20 ) #6 A3 1 20 ) # A A4 1 20 ) # NE 1 20 > # C A I 20 ) # CB I 
120 l 

Cl -,E*.SICn *lPh* I 2# 5 1 2 ) 

co-yon s#p-in#;-in#s:plt#sceror 

c x 4 n o y are the trac< and range coordinates on the imagery # 

C RESPECTIVELY. R AND S ARE The map POINTS AND SHOULD BE WITHIN 45 

C CE3RFES OF ThE I - A 0 £ R Y POINTS. Ni ANO N2 SHOULD LIE APPROXIMATELY 

C In T-E track DIRECTION *1Th ni TO The LEFT OF n2. 

CALL PlOTSI IB gF# 1000.7 I 

c plot tape on ias 

RE *C 1 1# 55 I AA1# a*2# aa3# AA4 
55 FOR- aT l 20*4 ) 

25 READ I 1# 1 JN#N1#n2/AA 

1 FORMAT I 315# 16A4 I 
IF i n I 20# 20# 30 

30 CONTINUE 

c scale of plot and scale of errors 
REaOI 1#3ISCPlT#SCER0R 
DO 2 I • 1 # N 

READ! 1#211 I P I I I # 3 ( I l 
211 FORyaTi 3F10.2I 

2 REaDI 1.3451X1 I )# Yl I I. ALPHA! 1# I I# ALPHA (2# 1 ) 

C 

c 

C C-AN3E format CARD 

• 3 rOR-ATiario.31 
3*5 FOR-*T(2rio.3,A4,A*) 

C RO T *TION of r*“ SECTION 

AaDATAN((TIN2)>YINX)l/IXIN2)*XINllM 

BaDATANllClN2l*0IM))/(PIN2)>P(Nl))l 

a«a»3 

1003 A’-Ol*** 

ACEa-57. 29578*4 
DO 4 I ■ 1 # N 
B*P I I ) 

P( I )a»( l ) a0C0SI •* 1431 I laDSlN! -A ) 

4 Cl I »«0( l la(OCOSI-A) l-BalDSlNI-A) I 
CALL LSF(1#N#T#X#3#P#A#B) 

DO 5 I ■ l # N 

5 El I laAaYl 1 l»B-CI I I 

CALL LSF(1#n#P#y#E#3#A#BI 

a»oaT an l a l 

*0E0aA0EG457 .29578*4 
"R I TE I 3# 800 I AOEO 

BOO FORHATllX# 'ANGULAR ROTATION *'#F10»4# • DEGREES ' I 

00 6 !• 1»N . 

BaPlII 

P| I laP| I I • COCOS I -A> ) *0 I I I * I D8 IN I »A ) I 


p °osjM 




C- 3 




APPENDIX C 


0 EDA -2 089 


l 


1 9 ! 3 5 APR 03 / ' 75 


5a. 

6 

St I)«3t 1 )*DCOSt-A)-a»DSlNl-A> 

55 • 

C 

PLOT SCALING section 

56* 


p":n»io«eio 

57. 


S-IN-10.E10 

56 > 


00 200 I a 1 J N 

59* 


IF t P 1 1 ) *L T *PMI •« >PWl Nap I 1 ) 

60 • 

200 

IF(2( II .LT.QMNiailNAQi I | 

61* 


CALL ®lOT ( 0 • # 1 • # -3 1 

62 • 


C Awl LSF( 1 aSaYaX,3aP,AaB) 

63. 


Ala* 

6* • 

1005 

BY-8 

65* 


A l • 8 

66< 


CALL LSF( i/N.X/Y/P/3#A/B) 

67. 

C 

TE-pqRaRY SECTION to MEASURE IMAGERY DISTORTION 

66 • 


A2a ( 41*4 | /2. 

69* 


00 51 1 • 1 . N 

70. 


XI ( I )-A2*X( I | 

7l« 

51 

»i t I )aA2*YI 1 1 

72- 


A 3«0 « 

73* 


A a a 0 , 

7*. 


00 52 I • 1 A N 

75- 


A3**3*Xl 1 I >-PI I 1 

76- 

52 

aa.u.t i (11*5111 

77* 


A3-A3/N 

78. 


A A a A a / >, 

79* 


DC 53 I • l / N 

80* 


Xl ( I |a*l ( I 1 a A 3 

81 • 

53 

Y 1 C I 1 • Y 1 ( I 1 a A A 

82* 


Call mCP t X 1 « V l a P* 3a N a AQa Va VI $ VE / XE* YEa C l / N I 

83* 


<<•0 

8* . 


»R I T£ < 3» 56 1 44l 

85* 

56 

FCRA aT(Im1///20Xa20AA///I 

86 • 


CALL EPLOT IPa3.XEaY£,EXaEYaAA1a AA ) 

87. 


GO TO 5A 

88* 

57 

CONTINUE 

89* 


DO 70 I" 1 a N 

90 > 


Xl ( I 1 ■ A • X < l I 

91 • 

70 

Y1 ( I 1 • A 1 • Y ( I | 

92. 


*3»0 • 

93* 


A A a 0 « 

9a. 


DO 705 I «1 a N 

95* 


A3*43axl t I )-P( I ) 

96* 

705 

A AaA a*YI (11*5(11 

97. 


A3*A3/n 

98* 


A A • A A /N 

99* 


00 710 I • l A N 

IOC* 


Xl ( I ) — X 1 ( I >-A3 

101 . 

710 

Y 1 1 1 I ay 1 ( I )*AA 

102* 


CALL hCP(X1aY1aPaBaNaAQaVaV1/VE/XE/YEaC1/M> 

103* 


wRITEt 3»56) 442 

10a . 


C*ll EPLQT<PaQ»KE«YEaEXaEYaAA2aAA> 

105* 


00 TO 5* 

106> 

67 

CONTINUE 


AFTER CORRECT ALIGNMENT 


<7 


r 

I 

L. 


If 
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107* 


C-.0005 

108* 


0*0. 

109* 


CO 7 l»l/*« 

110* 


E ( I ) *A*X ( 1 I ♦B* 9 t I I 

111 • 

7 

0«0*D*9S( £ 1 1 ) 1 

112* 


AL*C 

113* 


c-3. 

1 1*« 


00 8 I-l/N 

115* 

8 

XI 1 I )*X( I l*AL*Tl 1 > 

Ufa* 


Call LSF( 1/N/X1/X/P/3/A/B) 

117. 


CO 9 I-l/N 

118* 

9 

p»r«2A3S(A*xi(i>*a-P(i)) 

119> 


XF(F*D 111/11/ 13 

120 • 

10 

A L ■ • A L 

121 • 


c—c 

122> 

11 

*L« AL*C 

123* 


00 12 I • 1 / N 

12a . 

12 

XI ( I )-Xt 1 I*AL»r« I 1 

12S> 

1006 

BX*B 

126* 


CALL LS*< 1/N/X1/T/P/3/A/B) 

127. 


3«0 • 

128. 


00 13 I-IZN 

129* 

13 

3«3*0*as( *«Xl (I l*B“PC I ) 1 

1 3C • 


I r <G-*’ll6/16/l* 

131 • 

16 

*■■3 

132* 


GO To 11 

133* 

1* 

AL-AL-C 

13*. 

C 

T£ -ooRxRr SECTION 2 

135* 


A * • 0 . 

136* 


A 3*0 . 

137* 


CO 71 J»l/N 

138* 


TUI) «A2*T 1 I ) 

139. 


XI 1 I )*A2*X ( I ) 

1*0* 

71 

XI ( 1 >«X1 ( I >*AL*U t I ) 

1*1" 


00 72 I *1/ N 

1*2* 


A 3* A3*X 1 ( 1 )*P( I ! 

1*3* 

72 


1**« 


A3-A3/N 

1*5. 


a*«a»/n 

1*6* 


00 73 I • 1 / N 

1*7. 


XI 1 I ) ■ X 1 ( I 1 • A 3 

1*8* 

73 

T1 I 1 )«T1 ( I »-A* 

1*9. 


C*LL HCPIXl/Ti/P/Q/Nz AQ/ V/V1/VE/XE/YE/C1/N) 

150. 


■9ITE13/56IAA3 

151* 


Call EPLOT (P/Q.XE/TEz E Xz E Yz A A 3/ A A ) 

152. 


GO TO 5* 

153. 

68 

CONTINUE 

15*. 


DO 7* I • 1 / N 

155. 


XI ( I )*X( I l*AL*T( I ) 

156. 


XI I I )«A*Xl 1 I ) 

157. 

7* 

Tl ( r l»Al*T( I | 

158. 


A3«0 • 

159. 


A*«0« 
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19! 35 APR 03# *75 


1 6C • 


CO 713 I ■ 1 « N 


161 • 


*3a*3*xl t I l«P( I ) 


1 62 • 


713 A4a***Yl ( I 1 >3 ( 1 1 


163 • 


ASaAS/N 


164 • 


A.iA./N 


165. 


CO 714 laliN 


166. 


XI 1 1 )*X1 1 I )*A3 


167. 


714 X l ( I ) • Y 1 ( I |a*4 


166. 


CALL “CPlXljrliP.O.N.AQ^Y^Vl/VEiXE/YEiCl/tlt 


169. 


»RITC< 3 # 5 6 ) A A 4 


170. 


CALL E °L 0 T 1 P# 0# XE # Y £ i EX# £Y # AAA# AA ) 


171 . 


30 TO 54 


172. 


69 CCsT1\lE 


173. 


a L« a L • < *57.296 1 


174. 

C 

SCALE CARDS 30 "ERE 


175. 


-RITE 1 3# 35 ) A A 


176. 


»RITE(3#15IA#A1#AL 


177. 


15 FCRyaT ( lox# 'T hE T R a 0 < DIRECTION SCALE FACTOR IS 1!'#F12.4#' # 

THE 

176* 


1 RA \3E DIRECTIOW/lOX# 'SCALE FACTOR IS 1!'#F12.4,*. TnE SHEAR 

IS • 

179. 


2# F t, . 2 # ' OEGREES ' /// ) 


160 * 


35 FORMAT ( 10X# 16A4 1 


181 • 


30 TC 25 


182. 


20 COnT I s V E 


163. 

C 

ST t\o*RO "RITE SECTION 


184* 


30 TO 75 


185* 


54 *<a« <♦ 1 


186. 


«R ITE ( 3# 58 I 


167. 


58 FOR-AT ( 28X# ' TA9lE I ■ RESIDUAL ERRORS ' //17X# 'POINT NUMBER 1 #10X# 'TR 

18b< 


1 a C< ERROR '# 1CX# ' RAN3E ERROR'//) 


189* 


► RITE(3«60MI#XEIII»YE(I )# ALPHA 1 1# I )# ALPHA I 2# I 1 # >1# N ) 


1 9C • 


60 FOR* ATI IB*# 13# 1 6 X# F 1 3 • 2# 1 IX# Fi 0 • 2# 4X# A4# A4 ) 


191 • 


Sl"*»0 * 


192. 


SJ“TaO. 


193. 


CO 21 I ■ 1 # N 


194. 


S'J“aS , J"4XE I I ) axE I I 1 


195* 


21 SL'“YaSj-Y4YEt 1 )*TEI I ) 


196* 


SU-YaSLMY/N 


197. 


SU-aSUM/N 


198> 


SJIaSORT ( S um 1 


199. 


SJ*YaS3RT ( SuHY 1 


200. 


-RITE 1 3#22)SuH#Sj*iy 


201 . 


22 F0R“ATI5X#’X VARIANCE a * # Fl 0 » 4# 1 V VARIANCE •*/FlO»4> 


202* 


30 T0(S7#67#6B#69#7S)#KK 


203. 


75 CONTINUE 


204. 


CALL PlOTI 12«0<0.#999) 


205 • 


CALL EXIT 


206. 


END 
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l • 
2* 
3 • 
* • 
5 « 
6 • 
?. 
8 • 
9 * 
1C. 
11 . 
12* 
13 . 
1*. 
15 * 
16 * 
17 . 
18 * 

19 . 

20. 


19135 APR 03. 1 75 
SwSROWTJnE LSPl J* n*X* Y*P*Q*A*B> 

RE*L »8 X (1 l,Y( 1 I.Pl 1 l*Q(l I aS*R>T*U*A*S 
C T-IS SE C T I OS DOES A LEAST SQUARES PIT OP h POINTS* STARTING FROM 

C T-E JTh oojmT IS X ASO P 

A • 0 . 

s«c * 

R«C * 

e«o« 

T ■ C • 0 . 

L-0.0 

DO 1 I»J*K 
S-S*X(II 
R«R*P ( I I 
T»T»X 1 1 >•* 1 1 ) 

1 U-U*X(!>«P(I> 

A«(U«< (S«RI/Pt l/(T-( (S*S)/«) ) 

B» ( R-A*S ) /H 

RETURN 

END 


pa „ 


C-7 


APPENDIX C 


GERA-2089 


l • 
2* 
3* 

4 • 

5* 
6 • 
7 « 
8* 
9* 
1C' 
11 • 
12 ' 
13* 

14. 
15* 
16' 
17. 
16* 
19 > 
20 • 
21. 
22> 
23. 
2*. 


19:35 APR 03. • 75 

SwflBC'wTlNt E P L J T ( p . 3. X£ . YE.EX.EY. AAA. A A I 

CCjeLE p RtC IS I Os PI 1 t .,.( 1 I.XEtl). YE II I.Extl >.EYI 1 ) 

Cl -EsS I ON AAA ( 1 ) . A4 | 1 I 

CO--CS n. p ~In. ids. SC p lT.SCERQR 
C*lL Sy-bCl(C..0.4, .14, aaa.90*.80» 

X“AX*0. 

Call PlQTI 0.5. J..-3I 

OC 1 1-l.N 

x p aGE*(P<Ii*PhIM/SC b lT 

YPaGE* ( 0 ( I l*CM I M /SC B LT 

JE|XPA3E'OT.XflAX)X-AX»xPA0E 

CALL SYmbCli XPa 3E, Y B A3E, *07, 3. 0.,-t > 

x p a3E*x p aSE*x£I I I /SCEROR 

IP(xpa 3 E.qt.xiaxi xr ax*xpage 
yPaOE*YPa3E*Y£ ( I I /SCEROR 
1 CALL SymBOlI XP43E. Y b A3£, .07, 1,0.,*2> 
xPaQ£>xPaQE«.5 
ILIyi«Xpaq£ 
xPA0E“tLI f '**5 

CALL symbol (8* 000.1 .A, .140 a AA.90..6*) 

x p a0E*xpa0E*12. 

call plottage, o.,»ai 

RETURN 

END 


1 • 
2* 

3 • 

4 . 

5* 
6 • 
7. 
8 • 
9* 
10 * 
11> 
12 * 

13. 

14. 

15. 
16> 


19535 APR 03. '75 

SUBRCuTInE HCPIXl.Yi,P,O.N.AQ.V.Vi»VE/XE.YE/Cl.M) 

REAL >8 Xl i 1 ). Y1 I 1 ).»( 1 I/O I 1 ),A3I 1 ), VEI l >/ V/Vl.XEf 1 ). YE 1 1 ) 

c this SwSrcjtine com»jteb the mistoorah or the lenoth deferences between 

C T-£; Radar and -ap POINTS. AS well AS ThE ABSOLUTE ERRORS AND CEPICD 

CO 3 1*1.20 

3 aQI I l*c. 

00 4 I • 1 , N 

XE I I ) ■ X 1 ( I >*P< I I 

4 YE I 1 > *Y 1 I I I *0 ( I I 
CO 5 I-l.N 

ASI 1 )«AOI 1 )*DABSI XEIJ) » 

5 A 3 t 2 >»A 0 l 2 »*DABS( YEI I I I 
A< 3<1 loAQI 1 l/N 

A 3 I 2 )*AOl 2 ) /N 

RETURN 

END 


pa °z is 

1 POOB Q VAUTr 
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IMAGE DISTORTION CORRECTION 


D-l 






APPENDIX D 
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IMAGE DISTORTION CORRECTION PROGRAM 



*nma>BiQ aic s blank, tmmm» 
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Card input to the Image Distortion Correction program. 


Card 1 ITN 

(k values) ITO 

NAZ 
NRG 


Tape unit i nput • 
Tape uni t output . 
Number of azimuth 
Number of range e 


samples per record, 
ements . 


Card 2 SAZ = Azimuth scale. 

(3 values) SRG = Range scale. 

SKEW = Skew or shear angle in degrees. 
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c MAIN LINE PART Two INTENSITY T8 DENSITY/ 

C S * E w REMOVE POUTINE# SCALE IN AZIMUTH/ AND WEIGHT IN RANGE 

C 

DIMENSION ID { 21C0 > 

DIMENSION KIN (20) 

DIMENSION JD ( 2 1 CO ) 

C 

PI a 3 « 1 M 5926535897 
200 CONTINUE 


RF. AD 

99/ NN 



PRINT 

101 



PRINT 

102 



PRINT 

9 a / Nr-’ 



ITN 

TAPE U'-IT 

INPUT 


IT? 

TAPE UMT 

OUTPUT 

NAZ 

\U M -PER RE 

AZIMUTH SAMPLES per record 

NRG 

number er 

RANGE 

elements 

READ 

103/ ITN/ 

ITO/ 

NAZ/ NRG 

PRINT 

100/ ITN 

/ ITO# 

NAZ/ NRG 

3 AZ 

AZ I M(JTR 

SCALE 


SRG 

Range SCALE 


S u 't A 

SKEW OR 

SMEAR 

ANGLE 


READ 100/ SAZ/ SRG/ SKEW 
PRINT io?/ S A 7 / SRG/ SKEW 
REMIND I TP . 

PEW INC ITM 
RATIO ■ SAZ/SRG 

A l N ■ N A Z 
A ZN ■ AZ N *RATI 6 
RGN * NRG 

NZ ■ ABSIFG'' *SIN(3KEW*PI/l80. ) ) 

NaZP « AZN 

NAZP » NAZP+NZ 

BP«nZ 

A A a \RG 
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AA:.^P/\A 

IF (3KEW)220#22C#221 

220 COM I VJE 
5 13 * l • 

A A a - A A 

221 C^NTINUF 

call intcfn ud#-i> 

N2 = NZ-t-2 

C e h’F ON EACH SIDE ALWAYS 

DP 2C1 I RG « 1* NRG 
REAP TAPE ITN, ( ID(K),K-1,NAZ) 
CALL INTDEN{ ID#N'AZ) 

X I R 0 « I R 0 

I l* a A*XIRG*PB+.5 
IZ » PAX( !Z#1 ) 

JZ«nZ-1Z 

K« 1 

D9 222 I I » 1* I Z 

JD(<) =0 

K«<4 1 

222 CONTINUE 

DO 222 1 I • 1 # N A. 2 
wtMK)»IP(II) 

N = K. 4 1 

223 Continue 

D° 220 H«1*JZ 
JO ( < ) *0 

K * K 4 1 

220 C°NT I NUE 
N * N A Z ♦ N 7. 

CALL FILL! JO# I D# N# R AT 1 0 ) 

CALL RGv.GT <ID#IRG#N) 

WRITE TAPE IT9, (ID(K)#K«1/N) 

N « N I N ( 130# N ) 

PRINT 106# ( ID(K)#K«1#N) 

201 CONTINUE 

END FILE IT° 

REWIND 1TB 
PRINT 130 


sssias? 
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G r - TQ 20C 

98 F n R^ A T (6X2CAA//) 

99 rn^AT (PC A A ) 

100 F r R^AT ( ? F 1 0 • A ) 

101 FORMAT (6Xl4MFk'0r-RAM 9-75# 20X30HG00DYE AR AEROSPACE CORPORATION// ) 

10 2 FORMAT (feX33HSKFW CORRECT 1 0N AND SCALE PROGRAM# / ) 

103 F"RNAT( 1615) 

1 0 A Fr.RVAT{6XlOHl\PL/T TAPE/ 1 5# /6X1 IhOUTPUT TAPE/IW 

1 6X?8HNiJMp£P OF AZIMUTH SAMPLES I S# 15/ /6X27 MNUMbER OF RANQE ELEMEN 
2TS IS/ 16//) 

105 FORMAT (6X8HAZ SCALE/ F 1 2 • 1 / 5X8HRG SCALE Fl2« 1# 5X10MSKEW ANGLE Fl2*5 

1 //) 

106 F^RvaT (1X/13CI1) 

13A FORMAT (1P1) 

END 





noononooo 


A I'l'KNDIX 1 ) 


<TKI<A-2 08<) 


SUBROUTINE till <1D»JD#N,P) 

in IS THE INPUT ARRAY 
JD IS The OUTPUT ARRAY 

N IS THE NU^PE® OF P0 1 NTS IN THE INPUT ARRAY 
R IS ThE INCREMENT RAT 1 0 

\ WILL PE RETURNED AS THE NUMBER 0F 0UTPUT P0 1 NTS 


D I h r N s I m \ 10 ( ?1C0 ) 

DIMENSION JO (2100) 

I * 1 

X e 1 • 

X N = \ 

o * 2 

200 CONTINUE 

IE ( X-XN) 201# 201# 202 

201 CONTINUE 

X J 3 J 

IE(v-xu) 2C3# 2c3# PC* 

204 CONTINUE 
J * U* 1 
G'» TO 201 

203 CONTINUE 

0 v » ID(U)-ID( J-l ) 

Y « IC(J) 

Y - Y-DYMXJ-X) 

IE(Y) 2o5/ 205# 206 

205 CONTINUE 

JO ID » Y-.5 
GO TO 207 

206 CONTINUE 

JO ( I ) « Y^.5 

207 CONTINUE 
X ■ X*R 

I - 1*1 

GO TO 2oo 

202 CONTINUE 
N ■ 1*1 

return 

END 


(7 


Li 
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SUBROUTINE INTCf.N UD>Mi 
C INTENSITY T 0 DENSITY 

C 

D I K-F^S I fiN ID(?1C0> 

DIMENSION jr<P56) 

C 

C IF SCALE IS NOT KN0WN USE M.-N F0K FIRST CALL 

I F ( v ) 200/ ?1C# 210 

200 CO NTIKjr 

N * APS (K) 

K * ID ( 1 ) 

D*' ?C1 I * ?# N 
* * MX(K# 10(1)) 

201 CONTINUE 

<]$ - IMtNFITY SCAL0R USE 0 IF THE SCALE IS N0T KN0WN 

« density scalor use 128 for d»2 

RE AC 103# <IS# W'DS 
IF(KIS) 202# 202# 203 

202 CONTINUE 

c 

<is « <«6<* 

203 C(*K'TI\UE 

POINT 102# K I S# kDS 
al5F «= 1 • / AL0G ( 255 • ) 
xs * KDS 

DO ?CA I - 1# 25 6 
X * I 

D * ALGF«ALPG(X) 
w ir > ( I ) * D«XS*.5 

C ) I (D J #1 #301 TnIRP 

w0( l ) ■ 8 
cOA CONTINUE 
return 
c 

210 CeNTIN'JE 

D? 211 I • 1# r 
K * ID( I )/KTS 


K * NIN('<#2 c: ‘fe) 
K » Y A X ( 1 # K ) 

1PU ) « wP(K ) 

211 CONTINUE 
RETURN 


C 


10? F0RYAT ( 6X5HK I S 
103 F 0 R v A T ( 8 I 1 0 ) 
END 




• I 


IR#10X5WKDS .#110#/) 
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A IM'KNDIX I ) 


nKKA- 208 ‘) 


SUBROUTINE RGw'GT (ID,M,N) 

range weighting for vertical antenna pattern correction 

m a RANGE BIN 

\EG N Ff.tR INITIALIZATION IF REQUIRED 
DIMENSION ID ( 21 CO > 

IF (M ) 2pC, 201, 201 

200 CONTINUE 
RETURN 

201 CONTINUE 

return 

D» 202 I « N 
PUT weighting here 

202 CONTINUE 
RETURN 
END 
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